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ABSTRACT 

We study the propagation of uniformly translating fronts into a linearly unstable 
state, both analytically and numerically. We introduce a perturbative renormal- 
ization group (RG) approach to compute the change in the propagation speed 
when the fronts are perturbed by structural modification of their governing equa- 
tions. This approach is successful when the fronts are structurally stable, and 
allows us to select uniquely the (numerical) experimentally observable propaga- 
tion speed. For convenience and completeness, the structural stability argument 
is also briefly described. We point out that the solvability condition widely used 
in studying dynamics of nonequilibrium systems is equivalent to the assumption 
of physical renormalizability. We also implement a variational principle, due to 
Hadeler and Rothe, which provides a very good upper bound and, in some cases, 
even exact results on the propagation speeds, and which identifies the transition 
from 'linear'- to 'nonlinear-marginal-stability' as parameters in the governing 
equation are varied. 

Pacs numbers: 03.40.Kf, 68.10.Gw, 47.20.K 



1. Introduction 

The steady state equation for a travelling wave propagating into an unstable 
state does not uniquely determine the wave speed. As an alternative to directly 
solving the initial- value problem, velocity selection principles have been sought, 
which identify a priori the wave speed that would be dynamically selected from 
those permitted by the steady state equation. Such velocity selection and/or 
wave-number selection problems 1 14 arise in a wide variety of nonequilibrium 
systems that exhibit propagation of well-developed patterns and fronts into ini- 
tially unstable and homogeneous states. Examples occur in such diverse fields as 
population dynamics 15 and pulse propagation 16 in nerves in biophysical systems, 
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Taylor-Couette flows and Rayleigh-Benard convection in hydrodynamic sys- 
tems, crystal growth, 14 ' 17 models of solidification and aggregation, 18 as well as 
travelling waves in reaction-diffusion systems. 19 In the present paper, we consider 
only those systems for which the pattern emerging behind the propagating fronts 
is homogeneous, or can be described by an envelope curve, and concentrate on 
predicting the propagation velocities of uniformly translating fronts. 

1 2 20 

There already exist several proposed criteria ' ' for the dynamical velocity 
selection mechanism: a minimum speed rule, marginal stability and structural 
stability. Amongst them, the minimum speed principle has in certain cases a 
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rigorous basis in the Aronson- Weinberger theorem. According to the marginal 
stability hypothesis, 1 ' 2 ' 6,7 ' 8 for most sufficiently localized initial conditions, the 
propagation velocity of well-developed fronts generically approaches the marginal 
stability point which apparently coincides with the minimal velocity out of a 
family of stable propagating fronts. The marginal-stability point can be deter- 
mined explicitly from the linearized leading edge approximation, in which only 
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the linearized equation of motion is studied near the front. In the literature, 
this is sometimes referred to as the linear-marginal-stability case or the pulled 
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case. However, as emphasized by Ben- Jacob et al. and later by van Saarloos, ' 
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on the basis of rigorous work by Aronson and Weinberger for a class of simple 
equations, there exist cases in which the linear-marginal-stability selection fails. 
This is often referred to as the pushed case or nonlinear-marginal-stability case. 
Thus the marginal stability argument apparently provides a practical method to 
calculate the selected velocity analytically, but cannot tell when the method is 
reliable, since there is no general method known to distinguish a priori between 
pushed and pulled cases. This is an important limitation, because a given equa- 
tion may make a transition from pulled to pushed cases as a parameter in this 
equation is varied. 
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The most recent proposal is the structural stability hypothesis, in which 
the observable fronts are supposed to be stable against small changes in the gov- 
erning partial differential equation (PDE) itself. When the Aronson- Weinberger 
theorem is applicable, this hypothesis has the status of a theorem. The hypothe- 
sis applies (so far) without counterexamples to both the pushed and pulled cases, 
but until now it has not been able to yield an analytical means to obtain the 
selected velocity. Instead, a numerical method, based on the structural stability 
hypothesis, has been used. 

The principal purpose of this paper is to use the perturbative renormalization 
group (RG) theory to analyse the stability of uniformly translating fronts, and to 
show that this method can be used to calculate the change in front propagation 
velocity when the governing equation of a structurally stable front is perturbed. 
In addition, the RG allows us to predict the uniquely selected velocity by com- 
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bining the structural stability principle with it. As a by-product of our study of 
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fronts, we have investigated a variational principle due to Hadeler and Rothe 
for the phase trajectories of steady-state solutions of propagating fronts. We will 
see that the variational principle can give useful upperbounds on propagation 
speeds in both pulled and pushed cases, and is able to estimate the transition 
point between those regimes. This work, presented in section 6, is due solely to 
L.-Y. Chen. 

The applicability of RG to travelling fronts is predicated upon its recent 
successful application to the study of long-time global behavior of physical sys- 
terns described by nonlinear PDEs: porous medium equations, convection- 
diffusion transport with irreversible sorption, 29 the propagation of a turbulent 
burst, 30 ' 31 and linear continuum mechanics. 32 In these problems, the global solu- 
tions approach the similarity form u(x,t) = t~ a f(xt~P), where x is the spatial 
coordinate, t the time, and a and (3 are constants. In many cases - the so-called 
intermediate asymptotics of the second kind 33 - the exponents a and (3 can- 
not be determined by simple dimensional analysis. In fact, these exponents are 
analogues of the anomalous dimensions in field theory, and can be evaluated by 
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using renormalized perturbation theory. A mathematically rigorous formu- 
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lation has been given by Bricmont, Kupiainen and Lin; these methods have also 
been used to study front propagation in the Ginzburg-Landau equation. 35 ' 36 A 
detailed and pedagogical discussion of the application of RG to the asymptotics of 
partial differential equations, and its physical interpretation is given in the book 
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by Goldenfeld. We emphasise only that RG can be used to solve the above cate- 
gories of problems; standard techniques, such as multiple-scale analysis may also 
be used (perhaps with more difficulty), and in fact there is a general relation- 



ship between multiple- scale analysis and the RG, whose elements are discussed 
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elsewhere. 

A steady propagating wave solution has the form ip (x — ct) , where c is the 
propagation speed. If we introduce X and T as x = logX, t = logT, then 
ip (x,t) = \1/ (XT~ C ) for some That is, the solutions can be expressed in the 
form of similarity solutions. Because X and T must already be dimensionless, the 
velocity c can be viewed as an anomalous dimension, analogous to the exponents 
ct, (3 above in the intermediate asymptotics of the second kind. Thus it is a 
natural guess that there is an RG method to compute the velocity. However, 
there is an important difference between the similarity solution problems and 
the propagation problem. In the former, there is a unique asymptotic exponent, 
whereas in the latter the propagation speed is not unique. Thus, for the selection 
problem RG is not sufficient. We will see that it must be combined with some 
new physics — the structural stability hypothesis. 

The outline of this paper is as follows. In Section 2 a formal renormalized 
perturbation theory for the front propagation is formulated. In Section 3, techni- 
cal aspects of the renormalization method are discussed. The central issue is the 
nature of the eigenvalue spectrum of the linearised PDE; this is resolved using 
considerations of structural stability. In Section 4, we discuss the relationship be- 
tween solvability conditions and renormalizability. Section 5 contains a number 
of applications of the formalism, in which we compute the propagation velocity 
using perturbative RG. Section 6 presents results using the variational principle 
and a simple trial function. In section 7, we use the variational and RG methods 
to predict the transition between the pushed and pulled cases. We summarise 
and conclude in section 8. 
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2. Formal Theory of 
Perturbative Renormalization 

15 

In this paper we mostly concentrate on the so-called Fisher-KPP equation: 

where F is a continuous function with F(0) = F(l) =0. If F satisfies the 
condition F(ip) > for all ip G (0, 1), then there exists a stable traveling- wave 
solution interpolating between ip = 1 and ip = with propagation speed c for 
each value of c greater than or equal to some minimum value c*, where c* > 
c = 2a/-F'(0) (if F is differentiable at the origin). As suggested by the analogy 
between similarity solutions and travelling waves, there should be an RG method 
to compute the propagation speed. Here we show that when a front is structurally 
stable, we can devise a renormalized perturbation method to compute the change 
in the propagation speed due to a small perturbation in e.g. the 'reaction term' 
F. 

A. Perturbation Theory 

We consider formally a system described by the following abstract nonlinear 
(usually parabolic) equation: 

^= N W- (2-2) 

Let ipo(x — cot + xq) be a stable traveling front solution of (2.2) with speed Co, 
where xq is a free parameter indicating the translational symmetry of the problem 
explicitly; there is a one-parameter family of traveling wave solutions with a given 
speed. Let us add a small structural perturbation 5N to (2.2); for example in 



(2.1) we replace F with F + 5F. We assume that the operator norm of 5N on 
an appropriate domain is less than some small positive number e; for example 
for (2.1), the C°-norm of 5F is less than e: that is, \\SF\\ < e, where the norm 
\\SF\ \ = sup M \SF(u)\. Assume that in response the front solution is modified to 
ipo + 5tjj, where H^ll is °f order e. Linearizing (2.2) to order e in the moving 
frame with velocity cq, we find the following equation governing the first order 
correction 

^ = (c -^ + DN) ^ + 5iV{V> }, (2.3) 

where DN is the Frechet derivative of at ip = ipo, £ = x — cot + xq and the 
initial condition satisfies 6ip(£,to) = 0. 

To make our mathematical formalism as simple as possible, first, without loss 
of generality, we will consider the case for which the highest differential operator 
is the second order; however, the conclusions to be drawn below apply to other 
high-order nonlinear operators too. 

The formal solution to equation (2.3) reads 

t +oo 

<ty(£,f) = j dt ' j dt'G{^t-t\t')[5N{u Q m'), (2.4) 

to — oo 

where G is the Green function, which satisfies 



dG 
~dt 



- (c ^ + DN^j G = 5(t - t')5(Z - O (2.5) 
with G^0in|£ — £'| — > oo. The Green function may formally be written as 

oo 1 

GfotrfJ) = VSe^^T^^'^^^KCOeCt-O, (2-6) 

n=0 n 

where we use the summation symbol to represent both summation and integra- 



tion; is an appropriate weight function, {C^} are appropriate normalization 
constants, and 6 is the Heaviside function. The summation is over the spec- 
trum {cu n }, and {u n } are corresponding appropriately normalized (generalized) 
eigenfunctions obeying 

co-^r + DN) U n = ~U) n U n . (2.7) 



We do not need to solve completely this seemingly formidable equation. Differ- 
entiating equation (2.2) with respect to £, we have 

c°i + DN)( d 4f)=0, (2.8) 



dZ J\dC 

showing that, due to the translational invariance of the original equation, uo oc 
dtfjo/dt; satisfies (2.7). For the fronts in which we are interested, this function is 
well localized (square integrable), so that is an eigenvalue. We assume that the 
operator on the left hand side of (2.7) is dissipative, so that is the upper bound 
of its spectrum. 

Thus, only the zeroth eigenfunction •i/'o contributes to the secular term which 
is proportional to t — to, and the perturbed solution can be written as 

t) = MO - Sc(t ~to)^ + m r + 0(e 2 ). (2.9) 

Here 5c is given by 

f-™dtp(Ou mN{u m') 

and (6ip) r represents the non-secular terms, which are unimportant for the pur- 
pose of determining the modification of the velocity, leading only to finite 0(e) 
corrections to the profile of the front. 



One may immediately guess that this expresssion for Sc is the change in the 
front speed, but there are two problems with this identification. First, the naive 
perturbation theory is not controlled due to the secularity. A renormalization 
procedure, given below, will be used to remove the secularity. Second, both the 
numerator and the denominator of the expression for 5c may be divergent. This 
potential difficulty may be treated using the considerations of structural stability 
given in section 3. 



B. Renormalization and Renormalization Group 

For a certain class of sufficiently localized initial conditions, all the transient 
solutions to the unperturbed PDE of the form u(x,t) = u(x — x c (t, cq, xq), t), 
asymptotically converge to the same uniformly translating propagating front, 
u(x, t) = u(x — cot + xq) as t — > oo. After the perturbation has been switched on, 
xq is no longer a constant of motion of the perturbed system, and is therefore not 
observable at long times. From a measurement of xq at late times, the initial value 
of xq cannot be deduced; hence, it must be renormalized from the RG point of 
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view. As we saw in the introduction, if we make transformations x = log X, t = 
logT, xq = logAo, and u(x,t) = U{X,T), the travelling wave solution can be 
expressed in the form of a similarity solution U{X,T) = U (log AqX/T c °) = 
Ui {AqX/T 00 ) = U 2 (T/iAoX) 1 / 00 ^, for appropriate functions U, Ui and U 2 . It 
now becomes clear that the role of the constant of motion xq is the same as that 
played by the initial total amount of mass or energy Qq in the nonlinear diffusion 
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problems studied previously. Although it is possible to renormalize travelling 
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wave problems by regarding them as similarity solutions, here we will perform 
a new type of RG analysis directly on the travelling wave solution. 

9 



The divergence of the secular term (5u) s oc (t — to) as t — > oo, can be removed 
order by order in e by regarding to as a regularization parameter and introducing 
an additive renormalization constant Z = Z(to, fi, e) with fi an arbitrary time. 
Note that the limit t — to — > oo can be achieved in two equivalent ways: either 
keeping to fixed and letting t — > oo, or keeping t fixed and letting to — ¥ — oo. We 
will use the latter method in the following text, which corresponds to the con- 
ventional treatment of logarithmic divergences in the zero cutoff limit of the bare 
perturbation series for (e.g.) (p 4 field theory, t corresponds to logT in the analogy 
discussed in the preceeding section, so that t — to corresponds to log(T/Tb), where 
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log To corresponds to to- In the standard Gell- Mann-Low RG, we introduce an 
arbitrary length scale (or rather time scale in this context) L and split T/Tq into 
T/L-L/Tq; the term log(Tb/L) is absorbed into the (multiplicative) renormaliza- 
tion constant often denoted by Z . The corresponding procedure here is to split 
t — to into t — n — (to — A*) , where fi corresponds to log L in the standard approach. 
Since multiplications in the standard multiplicative renormalization procedure 
correspond to addition in our case, we renormalize xo(to) by Xq + Z(to, fi, e), and 
Taylor expand Z = Yl™=\ a n(to, fi)t n . The coefficients a n (n > 1) are determined 
order by order in e in such a way that all the secular divergences in ip(x,t) are 
cancelled out. In this way, the renormalized solution ip(x,t) remains finite even 
in the limit to — > — oo. Hence, to O(e), the solution ip(x,t) can be rewritten as 

ip(x,t) =V>o(£o) - 5c(t -t )^p- + 0(e), £ = x - c t + x 

=^ (^ + eai) -5c(t-t )^ + O(e), £ = x - c t + x$ (2.11) 

where O(e) refers to all finite terms of order e, regular in the limit to — > — oo. 
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By choosing a\e = Sc(/i — to), the secular divergence is removed. We obtain the 
renormalized solution 



t) = Vo (0 - Sc(t + 0(e), £ = x-c t + xf (//). (2.12) 



It is impossible that the actual solution ip(x,t) can depend on the arbitrary time 
scale [j,, because (as with L in the standard procedure) fx is not present in the 
original problem. This is expressed by the renormalization group equation 



dip 
dji 



= 0. (2.13) 

x,co,t,e 



Hence, to order e the RG equation yields, after equating \x with t, 



This has the form of an amplitude equation, an observation discussed further in 
section 4, and in ref. 23. Thus the speed of the renormalized wave is indeed 
c = Co + Sc, and the leading long-time asymptotic behavior, to 0(e), is 

ifj(x,t) ~ifj (x-ct + x )+O(e), (2.15) 

where x' Q is the new constant of motion for the perturbed system. 

C. Application to Fisher's Equation 

Now we consider the application of the above formalism to Fisher's equation. 
An important technical aspect of this discussion is deferred to section 3. We 
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consider only a perturbation SF to the nonlinear reaction term F of (2.1) so 
that (2.4) reads 



t +00 



S^ ,t) = e - c ^' 2 J dt' J d^G(^,t^',t')e^'/ 2 [8F{i; m'). (2.16) 



to —00 

Here, G is the Green function satisfying 

dG 



m -£G = 8(t-t')5(Z-a 



with G — > in |£ — £'| — > 00, and 



Formally, (7 reads 

= «o(o«S(0 +E e " <VB( '"° u »(0<(0, 

where £iio = 0, and £u n = X n u n . The summation symbol, which may imply 
appropriate integration, is over the spectrum other than the point spectrum {0}. 
Since the system is translationally symmetric, uq oc e c °£/ 2 V^(£). Due to the 
known stability of the propagating wave front, the operator C is dissipative, so 
that zero is the least upper bound of its spectrum. Hence, only uq contributes 
to the secular term in Sip. A similar argument as in the general case gives the 
explicit formula for (2.10): 

= G d^' (QSF{ip m (2 m 

ir^ co ^ 2 (o 
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3. Justification of Perturbative 
RG and Structural Stability 



Before proceeding further, we must now address the issue that the formula 
for Sc may not be finite, for example. The formal perturbation approach may 
fail if is not an isolated eigenvalue of the linearized operator. In this section, 
we will discuss when the restricted formula (2.17) is meaningful. To this end, we 
need results from our structural stability analysis, and so we begin with a brief 



Equation (2.1) can be classified into two cases: ambiguous and unambigu- 
ous. We say that the equation is unambiguous if it allows a unique propagating 
speed for its traveling wave solutions (here, we consider the waves traveling in the 
positive direction only), and is ambiguous otherwise. A necessary and sufficient 
condition for (2.1) to be unambiguous is that F does not have any isolated mini- 
mum at the origin. Notice that by a certain indefinitely small C°-perturbation 8F 
of F, we can convert an ambiguous equation into an unambiguous one. Hence, at 
most one propagating speed can be stable against this modification of the system. 
Because we are perturbing F, that is, the equation itself, we call such a perturba- 
tion a structural perturbation. Thus we may say that there is at most one struc- 
turally stable propagating speed for (2.1). It has been proved that the slowest 
propagating speed c* which allows stable (in the usual sense) propagating waves 
is a continuous functional of F so long as Q(F) = sup^ e ( ,i] F(ip)/ip is a contin- 
uous functional of F. This is satisfied, if 5F is C°-small and sup^w ,i] SF(ip)/ip 



is smaller than some positive number which converges to zero as \\5F\\ goes to 

zero. We call such perturbation a physically small (p-small) perturbation. Thus 
we may say that the wave with speed c* is structural stable against p-small 

13 



summary of this topic." 
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structural perturbations. 

We have conjectured that only structurally stable solutions of a model equa- 
tion correspond to the physically observable phenomena. This structural stability 
hypothesis is closely related to, but is not identical with the idea originally pro- 
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posed by Andronov and Pontrjagin for dynamical systems. The conjecture 
implies for (2.1) that c* is the selected propagation speed. It seems to be widely 
believed that c* is indeed the unique observable speed for (2.1) but a proof does 
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not exist for general F. However, our numerical studies have so far failed to un- 
cover counterexamples. The structural stability hypothesis is usually redundant, 
but has predictive power in cases where one's model of a physical phenomenon 
inadvertently includes unphysical features. This is the case for the over-idealised 
models (such as (2.1)) of propagation phenomena considered here. 

Prompted by our structural stability analysis, we have used the following 
numerical method to estimate the selected velocity. Apply the shooting method 
to the ordinary differential equation (ODE) governing the propagating wave front 
shape ip: 

where c is the propagating speed, and £ corresponds to x — ct, after modifying F 
by adding a p-small perturbation to convert the equation into the unambiguous 
class. We know that the unique speed is almost identical to the true selected 
velocity, since c* is continuously dependent on F. Therefore, there is only one c 
which allows the shooting method to give a solution. This c is the selected speed. 
However, in practice, the method may not be very accurate, because the true 
unique solution and other solutions which do not reach ip = for large £ may 
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not be easy to distinguish when the perturbations are small. 

There are several important consequences of our structural stability study. 
First, notice that even if 8F is p-small, — 8F need not be. Furthermore, it is easy 
to construct an indefinitely C°-small 8F which increases F'(0) indefinitely, but 
— 8F is p-small. An example of such 8F is a very tiny but very sharp spike well 
localized near the origin. This implies that for this SF, 8c must be indefinitely 
large, while 8c for — 8F is infinitesimal. Now (2.17) is a linear functional of SF, 
so that the sign change of 8F cannot cause such a drastic change. This clearly 
demonstrates that the formula (2.17) is not meaningful, if either 5F or — 5F is 
not p-small. More explicitly, sup u \8F{u)/u\ must vanish as \\8F\\ goes to zero. 
If 8F is differentiable at the origin, the condition is satisfied. 

The perturbation approach outlined above is legitimate only when is an 
isolated point spectrum of the operator C. Unfortunately, this is not the case 
for many reaction diffusion equations. For (2.1) the essential spectrum has the 
range (— oo, ??), where r\ = max{F'(l), F'(0)} — Cq/4, as is easily seen from Rota's 
theorem (a generalization of Weyl's theorem) 42 Hence, if cq = 2a/-F / (0), which 
is the important case for the pulled equations, then the eigenvalue is not iso- 
lated. Nevertheless, we can justify the formula or rather its augmented version 
as follows. 

Notice that we can always find a sequence {8fk} of piecewise differentiable 
perturbations such that 8fk converges to zero in C°-norm, but 8f' k (0) is always 
— 1. This is a sequence of p-small perturbations, so that c*(F + 8fk) converges 
to c*(F). However, the essential spectrum of the operators C is bounded from 
above by —1, so for these perturbed systems, is always isolated. Now consider 
the perturbed system with 8F. Instead of the system with F + 8F, we again 
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consider a sequence of systems F + 5F + Sfk- For this sequence, c* converges to 
c*(F + SF). Thus we may study the perturbation of the system with 5F + Sfk 
instead of the original system with SF to compute the change in velocity. Thus 
we may conclude that the ordinary perturbation theory can be used as done 
formally in the preceding section. 

For the pulled case, c* = 2a/-F / (0), so that both the numerator and the de- 
nominator of (2.17) diverge. This difficulty is removed if we consider the sequence 
mentioned above instead of the original system. In this case all the members of 
the sequence are unambiguous equations (pushed cases), so the decay rate of the 
eigenfunction belonging to is much faster than e 4 near the tip. That is, we 
have a natural regularizing factor for the members of the sequence. Taking the 
limit is thus equivalent to computing the limit with the aid of l'Hopital's rule: 



Thus with the aid of the structural stability consideration, we can give the 
correct form (3.2) for the change of velocity and a sufficient condition for its 
validity: ±5F must be p-small and the right and the left derivatives at the origin 
must exist. 



8c = — lim 



(3.2) 
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4. Amplitude Equations, Solvability and RG 

The general formula (2.10) for the modification of velocity due to the external 
perturbations can also be constructed from a solvability condition, analogous 
to that occuring in the dynamics of defects and dislocations in nonequilibrium 
patterns. 43 46 We start with the perturbed equation 



^ = N{^} + SN{^}, (4.1) 



which is assumed to have a steady-state solution of the propagating front 
ip(x,t) = ip(£), with £ = x — cot + xq. We further assume that the operator 
norm of 5N is of order e. Suppose that the perturbed velocity and solution can 
be written as c = c$ + 5c and ip = ipo + Sip, where cq and ipo correspond to the 
unperturbed velocity and solution. Then to order e we get 

- (DN + co Sip = 5c^ + 5N{iP }. (4.2) 

Because the linear operator on the left hand side of equation (4.2) has a zero 
eigenvalue corresponding to an eigenfunction uq oc dipo/d^, the condition for the 
existence of a nontrivial solution to (4.2) is that the right-hand side should be 
orthogonal to the null space of the operator: Vo = piO^o = PiO^uo/d^ where 
we have included the appropriate weight function p: 

5c (pVo, ^o) + (M>, SN{ip }) = 0. (4.3) 
We immediately recover the formula (2.10) for 5c. 
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In particular, let us now consider a case in which the change in the speed Sc 
results from the change of some parameter 7 in the nonlinear operator N(tfj,-f). 
If we assume that the traveling wave solution exists for some 7, and that we can 
write SN = (dN /d^Sj, and further assume that the order of the differentiation 
and the integration can be exchanged, then from (2.10) we obtain the formally 
exact result 



Unfortunately, we have not found a way to use this formula in cases where exact 
results are not already available. 

The reader may well wonder why we make a fuss about the RG in the preced- 
ing section, given that the formula for Sc is obtained trivially from the solvability 
condition of (4.2). Our point is to demonstrate that the solvability conditions 
(order by order) and the perturbative renormalization procedure are equivalent. 
The use of solvability conditions is reminiscent of pattern selection in dendritic 

14 

crystal growth phenomena, the derivation of amplitude equations, such as the 
time-dependent Ginzburg-Landau equation, or more generally the equations of 

47 

motion on the slow or center manifold (in the loose sense of the word). Indeed, 
(2.14) is the equation of motion governing the slow motion relative to the original 
unperturbed front. It is actually a general feature that amplitude equations and 
slow motion equations are RG equations. These general features of renormaliza- 

23 38 

tion group approach are discussed more thoroughly elsewhere. ' 



dc 
dj 



J^dte^jQdN/d^ 



(4.4) 
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5. Examples 

We now apply the perturbation approach to a variety of front propagation 
problems. 



A. Generalized Fisher's Equation. 

Consider the model: 

du d 2 u . , , 

Here, we identify F = u 2 (l — u) and SF = —jiu{l — u) and regard \x as a 
small perturbation parameter. The equation is unambiguous for all fi > —1/2. 
In particular, the unperturbed case is unambiguous and its unique propagation 
front shape is described by 

w = T77W7> £ = x-c t + xo, (5.2) 

where c* = c{ji = 0) = l/y/2 and xq is an arbitrary constant. From (2.17) 
we obtain 5c = —V2/j. Thus, the selected velocity of the perturbed equation, 
to O(fi), is c* = l/y/2 — \/2fi, which is the same as the exact result from the 
variational method we will describe below. Presumably higher order terms in fx 
will spoil this result. 

B. Van Saarloos' Fourth Order Equation. 

g 

Van Saarloos considered the following equation, 

d(j> d 2 (p d 4 (f) .^/...w, ,x /, ~\ 

m=d^-^ + b {b + <l)){1 -^ (5 - 3) 

where 7 < 1/12 and < b < 1. He numerically solved both this full time- 
dependent equation for ifj(x, t) and the corresponding ordinary differential equa- 
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tion for the steady propagating wave state in the moving frame with velocity 

g 

c for various values of b and 7. For 7 = 0.08 and b = 0.1, he observed that the 
velocity c ~ 2.715. 

Here we present analytical results obtained by treating the fourth-order 
term —'jd^cfi/dx 4 as a small perturbation. The unperturbed equation has the 
propagating-front solution 8 </>o(0 = (1 + e K ^) _1 , £ = x — cot + xq, where 
k = [cq — \Jcq- 4]/2, and c = 2 for 1/2 < 6 < 1; « = [c + y^j - 4]/2, 
and Co = V2b + l/y/2b, for < b < 1/2. The transition point from the pulled to 
the pushed case is b = b c = 1/2. It is worth mentioning that it is straightforward 
to obtain these results by applying the variational method presented in section 
6. From the formula (2.10) we obtain, for < b < 1/2, 

and 5c = -7 for 1/2 < b < 1. 

Note that the linear-marginal-stability velocity q switches from q = 2 to 
q 2 — 7 a: 2 VI — 7 as 6 crosses the value 1/2. We find that the 0(7) RG 
prediction agrees very well with the numerical calculation. For instance, for 

7 = 0.08, we obtain c = 2.696 + 0(7 2 ) for b = 0.1, which is close to the numerical 

g 

result c ps 2.715 by van Saarloos. 

C. Newman's Population Equation. 

The second example which we consider is a modified porous medium equation 

48 

which has been extensively discussed by Newman, in the contexts of population 
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genetics and combustion. The equation can be represented as 

du 1 d f ri du\ ^ . ^ n . 

ar=2ter&J +u(1 - !,) ' " £a (5 ' 5) 

For n = 0, it reduces to the Fisher-KPP equation and the propagation velocity 
is co = a/2- For n = 1, it has a unique travelling- wave solution 

ito(M) = (l - e? ) ©(-£)» £ = x-c t + x , c = ^, (5.6) 
where © is the Heavyside function. 

48 

Numerical results by Newman suggest that for any n > 0, solutions evolve 
asymptotically into travelling- wave solutions with a unique velocity c ~ (n + l) -1 , 
depending only on the value of n. 

To study this problem analytically, we write n = 1+5 and perform a perturba- 
tive RG calculation, regarding 5 as a small parameter. Expanding equation (5.5) 
in 5, we obtain, to 0(6), the perturbation SN{uq} = (6/2)d x (uod x uo\oguo). 
The weight function is = and the ground state is ipo(£) = duo/d^ = 
— e^, — oo < £ < 0. Using formula (2.10) we have 5c = —(13/48)5, and 
the perturbed velocity to 0(6) is c = 1/2 — (13/48)5. A consistent formula 
c = (2 + (13/12) 6 + 0(5 2 )y works remarkably well even for 6 as large as one. 
Setting 6 = 1, i.e. n = 2, we obtain c ~ 0.3243, which is in excellent agreement 

48 

with the numerical result c ~ 0.32 obtained by Newman. 

D. Pulled Case Revisited 

The renormalized perturbation theory result (2.10) can also be used to cal- 
culate heuristically the selected velocity of the unperturbed system by imposing 
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the structural stability principle. Within the perturbation theory, a necessary 
and sufficient condition that c* be the selected speed is that 5c(c*) must vanish 
as 5F vanishes. For concreteness, we consider the Fisher-KPP equation (2.1) 
with F = tp(l — ip). It is found that the change in the velocity 5c(c) is zero as 
\\5F\\ — > for all perturbations 5F, which are both p-small and differentiable 
at the origin, only for c = c* = 2; for c > c* there exist such perturbations for 
which 5c does not vanish as \\SF\\ — > 0. A simple example of the latter is the 
perturbation 5F = ip(l — ip) — (ip — A)(l — ip)Q(ip — A), where 9 is the step 
function, and we let A — > + . The RG calculation shows that 5c ~ \J c 2 — 4 as 
A — > + . Obviously, only for c = c* = 2, 5c goes to zero, but for c > c* = 2, 5c 
does not vanish. Therefore, only the wave with c = c* = 2 is structurally stable, 
and c = 2 is identified as the selected velocity of the unperturbed system. 

6. Variational Principles 

In this section, we digress briefly to report results of L.-Y. Chen using the 

24 

variational principle of Hadeler and Rothe. With relatively little effort, useful 
estimates for front speeds and transition points are obtained. 

We consider (2.1) in the variable u(x, £), where the reaction term F is contin- 
uously differentiable and satisfies the condition F(0) = F(u s ) = (where u = 
and u s are two steady-state solutions of the PDE); other generic conditions will 

24 

be imposed later. Hadeler and Rothe's principle is that 

c = inf sup + (6.1) 

P 0<u<u s P{U) 

where supremum means least upper bound etc., ' denotes differentiation with 

22 



respect to the argument, and the function p(u) = —d x u satisfies the conditions 
p(u) > 0, < u < u s ; p(0) = p(u a ) = 0, p'(0) > 0, p\u s ) < 0. (6.2) 

It was found empirically that if F is of the form F(u) = 11(11™ — u n )&(u), 
where $(tt) is some rational function which does not have zeros at u = or u s , 
then the following choice very often gives the correct c*: 

p(u) = u(u^-u n )i;(u,{K}), (6.3) 

where ip is a function positive for < u < u s and {k} is a set of variational 
parameters. Often one parameter k, which is the scaling factor of p, is sufficient. 
Even when the choice fails to give the exact c*, it gives a very good upper bound 
of c*. 

A. Fisher's Population Model. 

We consider the equation 

— = -^ + u(l-u)(l + vu n ), n > 1, -1 < v < +00, (6.4) 

where u = 0, 1 are two steady-state points. When n = 1, (2.7) reduces to the 
original Fisher's population model that Hadeler and Rothe discussed in Ref. 24. 
They indirectly obtained correct results for n = 1 by taking advantage of the 
known exact solution of the front profile although they did not choose a proper 
trial function. Here, we show that, for n > 1, the velocities and values of the 
parameter v at the transition from pushed to pulled case are the same as in the 
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n = 1 case, although we do not know how to obtain the exact solution of the 
front profile in any n > 1 case. In this sense, all n > 1 models belong to the same 
universality class, as far as the velocities and transition parameters are concerned. 
The numerical calculations which we have performed in several n > 1 cases by 
the shooting method are in excellent agreement with our analytical results and 
support the assertions above. 

We chose 

p(u) = ku(1 — w), (6.5) 

where k > is some adjustable parameter related to the decay rate of the front 
profile. If we define g(u) = p'{u) + F{u)/ p{u), then 



g ( u ) = K + - + -u n - 2ku. (6.6) 

K K 



Firstly, we consider the case v > and n > 1, and study the behavior of the 
function <j)(u) = uu u /k — 2ku with k > fixed. It has an extremum at the point 
uq = (2«: 2 /nz/) 1 / n_1 , which is determined by the condition 4>'(uq) = ui/Uq^ 1 / k — 
2k = 0, and further, since (j)"(uo) = n(n — 1)vuq~ 2 /k > 0, <fi has a minimum at 
u = uq. No matter whether uq lies in or out of < u < 1, the supreme bound or 
maximum is always reached at either of the two boundary sides, i.e. u = or 1. 
Thus, we have 



G(k) = max \g(u)} = k + - + max{0, - - 2k}. (6.7) 

0<u<l K K 



In the case of n = 1 and v > 0, we have g(u) = k + 1/k + {y j k — 2k)u, so we 
have the same result (6.7) as in the cases of n > 1. There are two separate cases 
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to consider below. If v/k — 2k > 0, i.e. k < \/v/2, then we have 



G(k) =k + - + (--2k) 

i + " K < 6 - 8 > 

= K. 



K 



Since G'(k) = — (1 + v) / k 2 — 1 < for any k, G(k) is a monotonically decreasing 
function of k, and the minimum is attained at k = (zz/2) 1 / 2 . Thus, in the interval 
v > 2, our estimate for the minimal velocity is 



c = min{G( K )} = ^. (6.9) 



If v/k — 2k < 0, i.e. k > then we have 

G(k) = k + - +0 > 2. (6.10) 

with equality only when « = 1. Also, we must satisfy the condition 1 > 
i.e. v < 2. Thus, the minimum velocity is 



c = min{G(K)} = 2. (6.11) 

K>0 



Secondly, we consider the case of — 1 < v < and n > 1. Because v/k — 2 k 

is always negative for any k > 0, we have the same results as equation (6.10) 
and (6.11). In summary, in all cases of n > 1, we have the following minimal 
velocities 

f 2, for -1 < v < 2 ; 

c =\ ^_ (6.12) 

1(2 + 1//^), forz/>2. 
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B. Other Examples 

A second more complicated example is the partial differential equation 

g 

For n = 4, it reduces to an equation considered by van Saarloos. By using the 
method of reduction of order of ODE, he obtained exact results for front velocities 
Co and the transition parameter value d c : 

J 2, for d < d c = 2/V3, 

co = < , (6.14) 

[ {-d + 2v / rf 2 T4/ v / 3), for d > d c . 

The nonlinear source term in the equation (6.13) can be written as 

F(u) = u{u n J 2 - u n ' 2 ){u n r /2 + u n > 2 ), (6.15) 

where uj 2 = [d + (d 2 + 4) 1 / 2 ]/2 and uj 2 = [-d + (d 2 + A) l / 2 }/2. It always 
has an unstable steady state u = and an absolutely stable state u = u s > 0. 
For simplicity, we only consider the propagation front connecting the above two 
steady states. 

In this case, we choose the trial function 

p(u) = Ku(u n J 2 - u n/2 ), k > 0, < u < u s , (6.16) 

but not p(u) = ku(u s — u) as in the former example. As a result, we obtain the 
complete result 



co = 



2, for d<d c = n/y/2(n + 2), 

[-nd +(n + 2)Vd 2 + A]/2^2(n + 2), for d > d c . 

(6.17) 



As a check, when n = 4, the result (6.14) by van Saarloos is recovered. 
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We remark that in these examples, only when g(u) reaches its maximum at 
u = u s > with k fixed is it possible for a transition to the pushed case to occur, 
while when g(u) reaches its maximum at u = it is the pulled case that occurs. 
We conjecture that this is typical. Based on this naive picture, we have also 
studied several other interesting cases: 

15 

(1) The Fisher-KPP equation, where F(u) = u — w. n , n > 2. If we choose 
p(u) = ku(1 — w. n_1 ), < u < 1, we obtain cq = 2. 

(2) Schlogl's second model for chemical reactions, 49 where F(u) = 7— /3u+3u 2 — 
u 3 = (l-u)(u-u 1 )(u-u 2 ),P = ^ + 2,0 < 7 < l, Ul = l + (l-7) 1 / 2 , u 2 = 
1 — (1 — 7) 1 / 2 . We choose p(u) = k(u — l)(u\ — u), 1 < u < u\, and find 



that the exact result cq = 3a/ (1 — 7)/2 is obtained. 

(3) A generalised version of Fisher's model, 24 where F(u) = u{l — u)(u — n), < 
u < 1/2. If p = ku(1 — u) is chosen, cq = l/y/2 — y/2/j, is obtained. This 
result is exactly the same as the perturbative RG result we obtained in 
Section 4. 

7. Transition from Pushed to Pulled Cases 

In this section, we consider how the transition point between pulled and 
pushed cases is changed when a p-small perturbation is present. 

Consider the following indicative example 

~qI = + u(u s -u)(u r + u + eu n ), n>0 (7.1) 
where u s = [d + y/d 2 + 4]/2, u r = [— d + V d 2 + 4]/2, d > and e is a perturbation 

g 

parameter. For e = the transition from the pulled to the pushed case occurs 
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as d is increased from zero at d = d c = l/y/2. We can apply the perturbation 
approach only for d > d c , because we need the explicit formula for the unper- 



turbed wave front: wq(0 = u s (l + Ae K ^) 1 , where k = [cq + Wc 2 , — 4]/2 and 



Co = [3V d 2 + 4 — d]/2y/2, and A is a constant of integration. For d < d c , there is 
no known formula for the unperturbed wavefront. 

The zeroth eigenfunction is w (£) oc — kuq{1 — uq/u s ) and the weight function 
p(0 = ( u s/ u o — 1) C °/ K . To determine the new transition point d c , we substitute 
into the formula (3.2) and obtain, to 0(e), 

_ u1+ 1 T(n + A)T(2-c Q /K) 

k; 3!r(n + 2-c /«) ' 1 J 

We also have c = 2 for n > and c = 2 a/1 + e-u s for n = 0. Setting c*(<i c ) = c(<i c ), 
we find that 



2 v / l + ew a (d c ) =c (d c ) +(5c(d c ). (7.3) 
For n = 0, we have 5c = \/2e and 

V2 + e 2 + 3e 

4 = s • (7.4) 



whilst for n = 1, we have 5c = e/2^j — 4 and 



4 - W 



For n = 2, we have 



5c = e/20 w s (c + 3^/c§-4)(co - a/c 2 , - 4) a/c 2 , - 4. (7.6) 
We can also apply the variational methods to determine the velocity c* and 
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transition point for any n > and d. We choose the trial function 

p(u) = ku(u s — u), < U < U s , K > 0. (7-7) 
For n = 0, we find d c = [y/2 + e 2 + 3e]/2 and 



{2^1 + w s e, for (i < (i c ; 

c = co(d) + y^e, for d > d c . 



(7.8) 



For n = 1, we have the result d c = (1 — e)/ a/2(1 + e) and 

2, for d < d 



C = 



^2/(1 + e)v / d2T4 + (e - 1)/^ + l)(d + Vd T T4)/2, for d > d c . 

(7.9) 

If we expand the above results to 0(e), we find that they coincide with the results 
from perturbative RG. We suspect that the results from variational methods are 
actually the exact ones for n = and 1. 

For n = 2, we obtain 

{2, for d < d C: 

^2/(l + e)(Vd T T4 + eu 2 s ) - y/(l + eu a )/2 u s , for d > d c , " ' U " 

where d c is the positive root of the equation eUg(d c ) +Ug(d c ) — 2 = 0. If we expand 
the velocities in (7.10) to 0(e), we observe that these results are slightly greater 
than those from perturbative RG, and yet give good upper bounds. In fact, for 
n > 2 and small e, perturbative RG gives more accurate results on velocities and 
transition points than the trial function we used. 



29 



8. Concluding Remarks 

In this paper, we have shown that for structurally stable fronts, a renormal- 
ization group method can be used to compute the change in the front speed 
when the governing equation is perturbed by a marginal operator; further, by 
combining the structural stability principle with RG, we are able to predict the 
uniquely selected front itself. Our results apply to both the pulled and pushed 
cases. We demonstrated that the solvability condition widely used in studying 
pattern selection in nonequilibrium systems is identical to the physical renor- 
malizability (observability) condition. We have also implemented a variational 
principle which gives very good upper bounds and sometimes exact results on 
front speeds, and which identifies the transition between the pulled and pushed 
cases. 

In future work, we will investigate the application of these methods for sys- 
tems where a spatial pattern forms behind the front. 
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